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Abstract 

Magnetic properties of 2D systems of magnetic nanoobjects (2D regular lattices of 
the magnetic nanoparticles or magnetic nanostripes) are considered. The analytical 
calculation of the hysteresis curve of the system with interaction between nanoob- 
jects is provided. It is shown that during the magnetization reversal system passes 
through a number of metastable states. The kinetic problem of the magnetization 
reversal was solved for three models. The following results have been obtained. 1) 
For ID system (T=0) with the long-range interaction with the energy proportional 
to r~P, the staircase-like shape of the magnetization curve has self-similar character. 
The nature of the steps is determined by interplay of the interparticle interaction 
and coercivity of the single nanoparticle. 2) The influence of the thermal fluctua- 
tions on the kinetic process was examined in the framework of the nearest-neighbor 
interaction model. The thermal fluctuations lead to the additional splitting of the 
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steps on the magnetization curve. 3) The magnetization curve for system with 
interaction and coercivity dispersion was calculated in mean field approximation. 
The simple method to experimentally distinguish the infiuence of interaction and 
coercivity dispersion on the magnetization curve is suggested. 

1 Introduction 

The properties of the magnetic nanoobjects and their systems are of current con- 
cern due to the appearance of technological possibilities of their fabrication and 
measurements. The reason of such interest is that such systems are ideal for study- 
ing collective effects and phase transitions and are also attractive as a media for 
high-density magnetic storage. Today there are the experimental data for the mag- 
netization curves of the 2D quadratic lattice of the Ni pillars with the perpendicular 
anisotropy 0, for the rectangular |] and square |]^, ^ lattices of the permal- 
loy nanoparticles, for the square lattice of CoCrPt particles with the perpendicular 
anisotropy for the system of the Fe |§] and permalloy ^ nanostripes, for the 2D 
systems of the chains of the Co particles [|10| , for anisotropic linear self-assembling 



mesoscopic Fe particle arrays |11, 12 



What is the main common feature of represented magnetic systems from the 
theoretical point of view? They consists of the magnetic coercive objects: mag- 
netic particles with the perpendicular single-particle anisotropy 0, 0]; particle 
chains, which have the effective anisotropy axis along the chain due to interparti- 
cle dipole-dipole interaction ||^, |TD|, [TT], |T^; magnetic nanostripes having the form 
anisotropy [§|, In the individual magnetic nanostripe the magnetization pro- 
cess take place by nucleation-propagation process |]T^ |T^. The propagation of the 



nucleus is very rapid. The numerical simulation demonstrates, that in ID nanopar- 



2 



tide chain the magnetization reversal proceeds through nucleation and the follow- 



ing domain wall propagation also []T5|. Both magnetic nanostripe and magnetic 
nanoparticle chain have two stable states within the external magnetic field with 
the magnetization directed along the stripe or the chain. A nanoparticle with the 
perpendicular anisotropy axis also has two stable states. The magnetization reversal 



in the particle has thermoactivated character [0. In the system of such magnetic 
nanoobjects there is the long-range interaction between them which have the mag- 
netostatic nature. In systems of the magnetic particles with perpendicular to the 
array plain anisotropy the interaction has the effective antiferromagnetic character. 
Its energy 

M(n)M(r,) 

= iTTs ' (1) 

M(rj) is a particle magnetic moment, Vij is interparticle distance. In the case of 
the system of the magnetic nanostripes the magnetostatic interaction are of long- 
range character also. It caused by magnetic charges appeared on the stripe edges in 
the magnetized state. The dependence of the interaction energy on the inter stripe 
distance is 

M(r.)M(r,) 1 1 ^ 

here M(rj) is stripe magnetic moment, is inter stripe distance, L is their length. 
For the neighboring stripes E ~ r~^, on the long distances E ~ r~^. 

Another system is 2D rectangular lattice of the magnetic nanoparticles with the 
single-particle anisotropy of the "easy plain" type. In this case particles form chains 
lying along the short side of the elementary rectangular. Due to anisotropy charac- 
ter of the dipole interaction the chain magnetization directed along the chain. The 



energy of the inter chain interaction involves the part (0) connected with the exis- 
tence of the magnetic charges on the chain edges and the part caused by discreteness 
nature of the chain W7\, ITBI]: 



Ediscr = ^ ^^^J^/ exp{-27rrij/a), (3) 

here M(rj) is magnetic moment of a chain, a is the interparticle (within chain) 
distance, Vij is the distance between the chains and L is their length. The relation 
Efiiscr/ E strip is proportional to the chain length. So for the sufficiently long chains the 



nearest-neighbor interaction plays a leading role The character of interaction 
is antiferromagnetic too. 

There were made some attempts to solve the problem of the magnetization pro- 



cess in the coercive system with the interaction by numerical methods [^, ^ 
It was found, that the magnetization curves look like staircase with the steps of dif- 
ferent widths. Besides there is the first report about the experimental observation 
of such steps ||T9 |. 



In our work we provide analytical investigation of the problem. In the first part 
we solve ID model of system of 2^ magnetic moments with coercivity and long- 
range interaction decaying proportionally l/r^ which corresponds to the system of 
magnetic nanostripes. We use the cyclic boundary conditions. It was shown that 
the magnetization curve consists of the series of steps corresponding to formation 
of superstructures. In the case N ^ oo this staircase-like curve has self-similar 
character (so-called "devil-staircase"). The method of the solution can be easily 
generalized on the case of 2D quadratic lattice of the magnetic nanoparticles with 
perpendicular anisotropy. In second part we solve the problem for ID system with 
niarest-neighbor antiferromagnetic interaction and single object coercivity at finite 
temperature. The influence of thermal fluctuation on the magnetization process 



described in first part is also discussed here. It is shown that the appearance of 
defects in the superstructures leads to splitting of the steps on the magnetization 
curve. The influence of the coercivity dispersion is also discussed. In the third part 
of the article we use the mean-field approximation to show, how one can distinguish 
interaction and coercivity dispersion in the system of magnetic nanoobjects. Very 
easy method to estimate the contribution of interaction and coercivity dispersion in 
the magnetic properties of the system is suggested. 

2 The model with the long range interaction. 
" Devil-staircase" 

Let us consider the ID system of long-range interacting coercive magnetic objects. 
Its appropriate example is system of finite length magnetic nanostripes or chains 
of magnetic nanoparticles. There is the effective antiferromagnetic interaction of 
magnetostatic nature between them. We consider its energy in the dimensionless 
form 

e = |, ^ |„ crfcO-n, (4) 
\k — n\P 

here ak — ±1 are interacting magnetic moments, n and k are the numbers of mag- 
netic moments positions, / is a dimensionless constant of the effective antiferromag- 
netic interaction (/ > 0). The nearest- neighbor object distance is equal to 1. Let the 
system is magnetized so, that all cr^ = —1. The magnetization reversal in the totally 
magnetised system onsets then the field at the object place exceeds its coercivity. 
In the case of the infinite chain it is 



N ^ 
k=l 



Here He is coercivity of the object and it is the same for all objects. The second 

term is a field originating at the place of one object by all the others. 

In a similar way it can be easily calculated, what reversal ends at the external field 

value 



H2 = H, + 2iap), (6) 

when the last object will be reversed. It is interesting to investigate magnetization 
curve in the region Hi < H < H2. As the system is one-dimensional, its ground 
state is disordered. Here we will consider the temperatures less than single object 
coercivity {kT < MHc), so the system can be in a number of metastable states. 
For example, if the interaction energy in the system is less than energy of coercivity 
there are 2" (n is a number of objects) metastable states at zero external field. 

If objects do not have coercivity or temperature of the system is rather high 
{kT » MHc) system will be in a ground state at any moment. If the external 
magnetic field is less than —2IC,{p) or larger then 21 ^{p) the ground state is totally 
magnetised state. The values of the magnetization of the ground states correspond- 
ing to external field at the interval —2IC,{p) < H < 2IC,{p) can be also found. This 
thermodynamical problem was solved in P2|. It was obtained that in this case 



magnetisation curve has steps and looks like self similar "devil-staircase". 

Here we solve the problem for the case when kT is less than any energy in the 
system. So system can be in metastable states and the problem of the magnetization 
reversal can not be solved by thermodynamical methods. We must solve the kinetic 
problem by correct choice of a consequence of metastable states the system passes 
through. The problem of such choice become easier in the case when the system 
consicts of absolutely even number (A^ = 2"') of objects and have cyclic boundary 
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conditions. 

So the magnetization reversal onsets then the external field exceeds Hi. Fluc- 
tuations are the reason that the place of the first reversed object can be choused 
arbitrary. Due to antiferromagnetic character of the interaction the additional in- 
crease of the external field is necessary to reverse second object. As the interaction 
decreases with the inter object distance, the second reversed object must be chosen 
at the largest distance from the first one. It is very easy to choose this place in the 
case of the system with cyclic boundary conditions (Fig. 1,1). 




^ l'^ x.x.x.x.x.x* 

I II 

Figure 1: The reversal process in the cyclic system of 8 magnetic moments and 
corresponding superstructures in the infinite system. Points mark not reversed mag- 
netic moments, crosses mark already reversed ones. The external field increases from 
a) to e). I) the reversal process in the cyclic system of 8 objects, II) the correspond- 
ing superstructures in the infinite system. 

In this case the magnetization process takes place through the sequential formation 
of different superstructures (Fig. 1,11) which are metastable states. Let us calculate 
the field values, then the definite superstructure appears {H_) and became unstable 
(i/+). At first we consider only superstructures with period m = 2^, k = 0,1,2... 
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and one reversed object per period (such in Fig. 1 b) c) and e)), k = corresponds 
to fully magnetized state. The field, when the superstructure with the period equal 
to m is formed (Fig. 2), is 



H(m) =H^ + 2I 



{2iap) - 21 



imP 



H,-2Ii{p)+AI 



(7) 



b b b b 

• X • • • X • K • X • • • X • • • • 

Figure 2: The process of the formation and following destruction of the superstruc- 
ture with the period m = 4. Points mark not reversed magnetic moments, crosses 
mark already reversed ones. H_ is the field when the magnetic moment in the point 
(a) becomes unstable and the formation of the superstructure takes place. if+ is 
the field when the magnetic moments in on of the points (b) become unstable. Due 
to fiuctuations one of them reverses first of all and then prevent others to reverse 
due to long-range antiferromagnetic interaction. 

Here the second term is the fields of already reversed objects. They prevent the 
chosen object to reverse. The term in brackets is the field of the other not reversed 
objects. They help chosen object to reverse. In a similar way, the field then the 
superstructure with the period equal to m become unstable (Fig. 2) is 



[m] 



He + 21 



(2^ - mp) 



mP 



- 2Ii{p) + 4/ 



(2/e(p) - 21 

(2^ - mv) 



{2P-l)i{p). 



mP 



(8) 



mP 



Here we use the relation 



OO 



E 



(9) 



The magnetization of the system is determined by the superstructure period and is 



, 1 ^ 2-m , , 

M= hm — y (Tfc = 10 

So there are the mast be steps on the magnetization curves corresponding to the 
stable superstructures, as the magnetization does not change while the magnetic 
field changes from H_ [rn) to Hj^ [rn) . Using (0,|[|lO]) we can rewrite the dependence 
of the step edges of the magnetization (Fig. 3) in the form: 

H^ = H,- 2Ii{p) + 4/e(p)(^^)^ (11) 
H^ = H,- 2iap) + 4/(2^ - m{p){^^r, (12) 
M 



1 


m=1 






o-_ 


m=2 








/-^'m=4 


-^- 


: m=8 



H 



Figure 3: The steps on the magnetization curve corresponding to the simple su- 
perstructures (m = 2^) with one reversed object per period. Magnetization of the 
corresponding step is M = (2 — m)/m = —1 + 1/2^^"-^ 

The steps corresponds to the reviewed superstructures do not cover all fields be- 
tween values Hi and H2 (Fig. 3). To analyze the magnetization behavior of the 
system while it transfers from one step to another we must take into account the 
formation of the more complex superstructures (Fig. Id). Let us consider, for ex- 
ample, the magnetization process between iJ_|_(m = 2) and H2 (Fig. 3), that is how 
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antiferromagnetic superstructure becomes fully magnetized. The superstructures 
with periods m2 = 2, 4. ..2' (/ = 1,2...) appearing with antiferromagnetic (mi = 2) 
ClS db background are represented at Fig.4. 



a) •« « • •«•••«•••«• 

mi = 8 mi = 4 

b) nXn n n nXn iO<» . lOC* ■ »X» 

mi = 2n^=8 mi = 2 m2=4 

Figure 4: a) are the simple superstructures (mi is a period). Points mark not 
reversed magnetic moments, crosses mark already reversed ones, b) are complex 
superstructures characterized by two numbers mi and m2. Small crosses mark 
antiferromagnetic background, big crosses mark reversed magnetic moments forming 
period of a superstructure. 



mi = 2,m2 = 2 is totally magnetized state, mi = 2,m2 = oo is antiferromagnetic 
superstructure. In this case M = 2/m2. The expressions for the magnetic fields 
H_{m2) and if+(m2) differ from the case of the simple superstructures as we must 
take into account the field of the antiferromagnetic background affecting on the 
objects. This additional field is 



,(2^ -l)e(p) _,^ (2^-l) e(p) 
ml 2P 



H = 2/^^^-^ = 2/^^^,^ (13) 



We must take this field into account twice as formerly it was directed along the 
external field but now it is directed against it. So we have 



i/_(m2) = H,- 2Ii{p) + mP){^r + 4/^^^^ (14) 
i7+(m2) = H,- 2iap) + 4/(2^ - l)e(p)(f )^ + 4/ ^^'~^/^^^ (15) 



Evidently (n,T3|) are exact similar to (|rT],|T^, but this curves starts at the point 
M = 0, H = H+{mi = 2) = H^- 2/^(p) + 4/ ^^""^^^^ which is the right edge 
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of the step corresponding to the antiferromagnetic structure, instead of M = — 1, 
H = Hi = He — 21^{p). In a similar way all other steps on the magnetization curve 
can be obtained. Each step corresponding to the superstructure is the base of a 
series of the steps, corresponding to the more complex superstructures with the first 
one as a background. The dependencies of if_ and -ff+ in all cases are similar, and 
H± ~ M^. So the picture become self-similar (Fig. 5). 

M 



0- 



m 



=2 m2=4 m3=8 



■ mi=2 m2=8 



mi -2 



H 







=4 m2=8 






rrii =4 








rrii =8 











Figure 5: The self-similar "devil staircase". 

Complex superstructures are characterized by a series of numbers mi, 7722,7/13..., 
where mi = 2'^ and rrii < rrii^i. Maximal rrimax is the period of superstructure, 
mi, m2...777maa;-i characterizes its background. Narrow steps are between wider ones. 
Let as find the common width of all steps. The step width can be easily calculated. 



2^ — 2 

AH = H+{m) - H_{ni) = 41— — ^{p) 

TTlmax 



(16) 



the number of the steps with equal width depends on m-max and = TUmaxI^. The 
common width is 
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oo oo op n oo op o 

^H= Y: N{n)AH{n)= ^ 2^/^^e(p) = 4/e(p) E 2'^'^^ 
n is even n is even ^=1 

(17) 

= 2/e(p)(2^ - 2) E(^)' = 2/e(p)(2^ - 2)^3117^ ^ ^^^^^^ 
So the sum of all steps widths is exactly the same as the width of the whole incli- 
nation of the magnetization curve, which is equal to 4J^(p) according to 

The difference in the step width is connected with the decaying character of the 
long-range interaction. The wider ones are conditioned by interaction of more near 
magnetic moments; the narrow ones are caused by the interaction of more distant 
magnetic moments. 

The magnetization values corresponding to the steps are 

4m — 1 , , 

M= 1, (18) 

where m = 1, 2...mmax/2. 

It is interesting, that the exponent in ([Tl|) and ([121) is equal to the power index 
in the dependence of the interaction on the inter object distance (H). So it is pos- 
sible to find the power index for the interaction by the experimental measuring of 
magnetization curves. 

There are some reasons of the distortion of the represented ideal picture. At first 
it is thermal fluctuations. As for the long distances {kT > Eint{r), k is a Boltzman 
constant) it leads to fuzzifying and disappearing of narrow steps which are condi- 
tioned by interaction of distant magnetic moments. Besides thermal fluctuations 
can lead to appearance of the defects in superstructures. This will be discussed in 
the next section. The second reason is the bounds of the real system. They can 
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play a significant role as the interaction has a long-range character. Nevertheless if 
the dimension of the system is larger than tq {kT = i?j„t(ro)), the influence of the 
bounds will be neglected by thermal fluctuations. Lastly the dispersion of the ob- 
jects coercivity can dramatically change the magnetization curve. Such self-similar 
behavior can be observed only in the system with small (less than interaction) co- 
ercivity dispersion. The method how to distinguish the influence of the interaction 
and coercivity dispersion in possible experiment is discussed in the last section. 

In spite of all deficiencies of the proposed model it help to understand the pe- 
culiarities of the magnetization process in the system, the nature of the steps on 



the magnetization curve ||T9|, ^ ^ and especially the fact that the difference in 
steps width is a sequence of the decaying character long range interaction. It also 
makes understandable the fact of the alternation of the narrow and wide steps on 
the magnetization curve [|T^, It is very probably that in the general case the 
magnetization curve has self-similar character too. The model can be easily general- 
ized on the case of the 2D square lattices of the magnetic nanoparticles. In this case 
one must provide the summation of the dipole sums for the superstructures on the 
square lattice. The superstructures must have square elementary cell in this case, 



so the dipole sums can be easy calculated |24 



3 The nearest-neighbor model. 
Thermal fluctuations 

Let us consider the magnetization process at finite temperature less than the coercive 
energy of a single magnetic moment kT < H^M but higher than interaction energy 
of distanced magnetic moments. In this case magnetic moments at distances larger 
than To {kT = E{rQ)) begin to reverse independently as their interaction is smaller 
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than the temperature. Nethertheless, as kT < HcM system can be in metastable 
states. The kinetics of magnetization process can be quahtatively described as fol- 
lows. When the external magnetic field exceeds the value of Hi the thermoactivated 
reversal of the individual magnetic moments begins. But already reversed magnetic 
moments prevents neighboring magnetic moments (lying at distances smaller than 
ro) to reverse due to effective antiferromagnetic interaction. More distant magnetic 
moments can reverse, as the interaction energy in this case is smaller than tempera- 
ture. Magnetization process has Poisson-like character and ends when the distance 
between reversed magnetic moments will be in the interval tq < r < 2rQ. The ad- 
dition external field is necessary to overcome antiferromagnetic interaction and to 
continue the process of the reverse. Evidently the thermal fluctuations will lead to 
distortion of the ideal picture described in the previous section. It is difficult to 
take thermal fluctuations in to account in the general case. Here we have solved 
the problem for the situation when the temperature is large than any interaction 
in the system with exception of the most powerful nearest-neighbor interaction. In 
this case the problem can be solved in the nearest-neighbor approximation [p = oo 
in (^). This model is also appropriate for the case of the planar system of the long 
chains of the magnetic nanoparticles when the main term in the interaction is inter- 
action of nearest neighbors (^). The form of the hysteresis in this case is represented 
in Fig. 6. The magnetization reversal starts at the field value H = He ~ 21, as the 
antiferromagnetic interaction helps the external field. 

The chains begin to reverse the magnetization due to thermal fluctuations. But if 
the one chain reverses, it prevents the neighboring chains to reverse as the effective 
field of the interaction is opposite to external field. Due to chaotic character of the 
sequence of the magnetization reversals of the elements of the system the defects 
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Figure 6: The hysteresis loop in the case of the nearest-neighbor interaction, a) 
in the case of T = the step corresponds to antiferromagnetically ordered super- 
structure, b) the dividing of step in two due to formation of defects in the case of 
thermal fluctuations. Mst is a step high. 

appearances is possible (Fig. 7), and the antiferromagnetic ordering with M = 
does not be achieved at this value of the external field. The additional external field 
[H = He) is necessary to reverse the defects. Then defects change their sign (Fig. 7). 

a a 
b b 

Figure 7: (a) are defects appearing due to Poisson-like character of the magnetiza- 
tion process, (b) are reversed defects. Points mark not reversed magnetic moments, 
crosses mark already reversed ones. 

The magnetization reversal ends within the field H = Hc + 2I when system becomes 
totally magnetised. So two steps appear on the each branch of the hysteresis loop. 
Their width is AH = 21 i.e. it depends on the interaction value. The magnetization 
value Mst corresponding to the step depends of defects concentration. It is the 
special problem to find this concentration. Let us consider the kinetics of the defects 
appearance in the system of N objects. Firstly all objects are magnetized against 
external magnetic field cr^ = —1. When the reversal process begins, the reversed 
objects begin to divide system into the regions of the not reversed objects (we will 
name them as "domains"). In process of magnetization the number of domains 
increases, their widths decrease. Let P„ is a number of domains consisted of n not 
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reversed objects then 



N-2 N-2 

M=T.Pn-T. ^Pn, (19) 
n=l n=l 

If a is the probabihty of the chain reverse per unit time, then 
dF ^ 

_^ = _ain-2)Pn + 2a ^ P^, n > 2. (20) 

k=n+2 

The first term is decrease of the domain number due to its dividing into smaller 
ones; the second term describes the appearance of new domains due to dividing of 
wider ones. As the domains consisted of one or two objects can not further divide, 
their number increases only. So 

dP ^ 

^ = 2a Y: Pk, n = l,2. (21) 

k=n+2 

Evidently a depends on temperature, but as we in the states stable at t — ^ cx), a 
does not affect final result. It can be easily checked, that 

d ^ 

Y.in + l)Pn = 0, (22) 

'^'^ n=l 

i.e. the whole number of objects in the system is constant. In the course of time 
{at » 1) only the domains with n = 1,2 remains. So P„(t oo) = for n > 2 
and M = —P2{t oo). Let us use the Laplas transformation and define 

POO 

Puis) = / Pn{t)e-''dt. (23) 
Jo 

Then, according to (|20D 
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where 



N 

Qn{s)= E Pkis). (25) 

k=n+2 

For s ^ 



AT N 2 

Qn{0) = E PdO) = E 7— ^Q'^(O)- (26) 

k=n+2 k=n+2 

In the recursive form this equation can be written as 

Qn(0) =Q„+i(0) + -g„+2(0) (27) 

n 

To find the magnetization value corresponding to the step it is necessary to calculate 
P2{t oo). Evidently 

P2(t oo) = limsp2(s) = 2aQ2(0) (28) 

s— >0 

If the maximum number of the objects in the system is A^, (^2(0) can be found from 
the (l27|) with the initial conditions Qn~2 = Qn-s = Pn{0), which are the sequence 
of (p^). In its turn pAr(O) = l/a{N — 2) according to (|20|) . The solution of the 
recursive equation was found numerically as 

Mst = 2a lim ~ 0.134 (29) 

N~*oo N 

As Q2{0,N) is proportional to a^^, the result does not depend on the value of a. 
So the defects formation during magnetization leads to appearance of two steps 
on the magnetization curve in the case of the nearest-neighbor interaction. The 
fiuctuational character of the magnetization process prevent chance to find system 
in the antiferromagnetically ordered ground (if H = 0) state. One can expect in the 
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case of the long range interaction thermal fluctuations will lead to similar splitting 
of the magnetization steps too. 

4 Dispersion of the coercivity 

Another reason of the distortion of the described magnetization process is the dis- 
persion of the coercivity in the system. It is necessary to distinguish the effects of 
the interaction and the coercivity dispersion especially in experiment. This problem 
is solved here in the mean field approximation. In the framework of this model the 
interaction is independent on the distance and e = I/N. It must be mentioned that 
this model is appropriate in the case of the very long stripes when L > Na {L is 
a stripe length, a is interstripe distance). The dependence of the magnetization on 
the external field is linear in this case as 

H^Hc + JM (30) 

He is coercivity, JM is the field of interaction in the mean field approximation. The 
main feature of the magnetization process in coercive system with the antiferromag- 
netic interaction is connected with its multistability. It means that if we change 
the sign of the external field changing the system does not change its magnetization 
immediately. At first it transits through the whole hysteresis loop from one branch 
to another (from point A to point B, Fig. 8), 

and then begins to change according to the new branch. In this case \Ha — Hb\ = 
2Hc. It must be noted that such behavior does not depend on the interaction man- 
ner. So different states of the system (characterized by the different magnetization) 
correspond to the same value of the external field. Such type of the multistability we 
will name interaction- type (I-type) , because the multistability can be in the system 
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Figure 8: The hysteresis curve in the mean field approximation, a) in the case 
of the system with interaction, b) in the case of the coercivity dispersion (uni- 
form distribution of the coercivity). c) system both with interaction and coercivity 
dispersion. 

of the non-interacting magnetic objects with different values of coercivity also. In 
this case magnetization reversal begins when field reaches the value Hi = Hcmin, 
when the reversal of the objects with the smallest coercivity starts. The reversal pro- 
cess is finished at the field value H2 = Hcmax, corresponded to the largest coercivity 
in the system. The hysteresis loop in this case the similar to the one for the system 
with interaction (as its branches can have similar inclination in the case AHc = 2/ 
and uniform distribution of He), but the transitions inside the loop differs. If one 
change the direction of the reversal process in the A point (Fig. 8) in this case, the 
magnetization does not change while the external field reaches the value of —Hi 
when the reversal of the objects with the smallest coercivity happens. This type 
of the multistable behavior we will name coercivity- type (C-type). The hysteresis 
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loop of the system both with the interaction and coercivity dispersion can be easily 
calculated too. It is represented on Fig.Sc. So one can distinguish interaction and 
coercivity dispersion by the behavior of the magnetization within hysteresis loop by 
analysis of 

the multistability type of the system. Evidently the self-similar behavior of the 
magnetization can be experimentally observed in the systems with small dispersion 
of coercivity, that is in the systems which demonstrate I-type of multistability. Let 
us note, that in the reviewed experimental works the type of the multistability is 
not examined in spite of the simplicity, from the on hand, and, significance, from 
the other hand, of such investigation. 

5 Conclusions 

By means of the simple models we have investigated the magnetization processes 
in the systems of the coercive magnetic objects with interaction. The reason of the 
formation of the steps on the magnetization curve is investigated. It is shown that 
the magnetization curve can have self-similar character. Its form is calculated in 
the case of the long-range interaction with E ~ The influence of the thermal 

fluctuations is analyzed in the framework of the nearest-neighbor approximation. It 
is shown that fluctuations lead to splitting of the steps on the magnetization curve. 
The effects of the interaction and coercivity dispersion on the hysteresis loop are 
examined by mean field approximation. Their difference is shown. It consists in the 
character of the magnetization dependence on the external magnetic field during 
the transition between the branches of the hysteresis loop. The understanding of 
this difference is very important to interpret of the experimental data. It allows to 
experimentally distinguish the influence of the coercivity dispersion and interaction 



20 



in the system on the magnetization processes. 
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